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Abstract 

(D 

' This paper is concerned with the dynamics and interactions of Q-balls in (1+1)- 

dimensions. The asymptotic force between well-separated Q-balls is calculated to show 
\ that Q-balls can be attractive or repulsive depending upon their relative internal phase. 

1/^ ' An integrable model with exact multi-Q-ball solutions is investigated and found to be 

On ■ of use in explaining the dynamics in non-integrable theories. In particular, it is demon- 

, strated that the dynamics of small Q-balls in a generic class of non-integrable models 

\ tends towards integrable dynamics as the charge decreases. Long-lived oscillations of 

' a single Q-ball can also be understood in terms of a deformation of an exact breather 

QO . solution in the integrable model. Finally, we show that any theory with Q-ball solu- 

tions has a dual description in which a stationary Q-ball is dual to a static kink, with 
an interchange of Noether and topological charges. 



1 



1 Introduction 



Q-balls are time-dependent non-topological solitons which carry Noether charge associated 
with a global U{1) symmetry of a nonlinear field theory [1]. Although they arise in a variety 
of theories perhaps the most important example is the minimal super symmetric standard 
model, where the Noether charge is associated with the U{1) symmetries of baryon and 
lepton number conservation and the relevant fields correspond to squark or slepton particles 
In this scenario Q-balls are condensates of a large number of squarks or sleptons 
and could play a role in baryogenesis through the Affleck-Dine mechanism [Ij. Potential 
interesting cosmological consequences include contributions to dark matter and isocurvature 
baryon fluctuations. 

In order to fully investigate the properties of Q-balls a crucial ingredient is to understand 
Q-ball dynamics and multi-Q-ball interactions. This is not a simple task as Q-balls have 
time-dependent internal phases and are solutions of nonlinear fleld theories. 

There is a vast literature on stationary Q-balls (see the recent work [19j and references 
therein) but only limited studies on dynamical aspects of Q-balls, mainly involving numerical 
simulations [3[|2]. These numerical investigations reveal complicated phenomena including 
phase-dependent forces, charge exchange and Q-ball flssion and fusion. These results are 
qualitatively similar for a range of theories and spacetime dimensions, including (1-1-1)- 
dimensions. Despite this fact, there is very little analytic understanding of Q-ball dynamics 
in any system. 

In this paper we address this issue by performing some of the flrst analytic studies of 
Q-ball dynamics. We study Q-balls in (l+l)-dimensions and compare our flndings with 
numerical results. We calculate the force between well-separated Q-balls, verifying that the 
force can be either attractive or repulsive depending upon the relative internal phase be- 
tween the Q-balls. We make extensive use of an integrable model for Q-ball dynamics to 
understand Q-ball interactions and perturbations in non-integrable theories. Furthermore, it 
is demonstrated that the dynamics of small Q-balls in a general class of non-integrable the- 
ories is captured by the integrable theory, with increasing accuracy as the charge decreases. 
Finally, we describe how any theory with Q-balls has a dual description in which a stationary 
Q-ball is dual to a static kink, with an interchange of Noether and topological charges. 

Part of our study involves a generalization of Q-balls to theories with a non-standard 
kinetic term. It is pointed out that a certain class of such theories is more conducive to 
Q-ball existence than standard theories, as most of the usual constraints on the potential are 
removed. This is true in all spacetime dimensions and therefore suggests new possibilities for 
Q-balls in systems previously not considered. Furthermore, non-standard kinetic terms are 
currently popular in cosmological theories and string theory, so Q-balls may have applications 
in these contexts. 
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2 Q-balls 



In this section we first review some material regarding Q-balls in theories that are typically 
used as the simplest examples possessing Q-ball solutions. We then turn to a discussion of 
modified theories with a non-standard kinetic term and explain why Q-balls are more generic 
in a class of theories of this type, since most of the usual restrictions imposed on the potential 
are no longer necessary. 



2.1 Q-balls in standard theories 

Q-balls arise in various types of field theories, but the simplest example consists of single 
complex scalar field with a global U{1) symmetry. Explicitly, consider the (1-1-1) -dimensional 
theory with Lagrangian density 

c = d,<pd>'^-um, (2.1) 

where, without loss of generality, we have chosen to set U (0) = 0. 

The Noether charge associated with the global U{1) symmetry is given by 

/oo _ _ 

(0(9*0 - 09*0) dx. (2.2) 
-oo 

A Q-ball positioned at the origin is a stationary solution of the form 

<P = e'^'f{x), (2.3) 

where > is the internal rotation frequency and f{x) is a real profile function satisfying 
the boundary conditions ^(0) =0 and / — as |x| — oo. 

Of course, given the stationary solution (12. 3p the position can be shifted by a spatial 
translation and the Q-ball can be given an arbitrary velocity up to the speed of light (which 
is one in our units) by applying a Lorentz transformation. 

Using the stationary ansatz (12. 3p the Noether charge takes the form 



/oo 
fdx. (2.4) 
-oo 



The classical field equation which follows from the Lagrangian density (12. ip is 

d,dy + ^ = 0, (2.5) 

which has a solution of the form (12.30 providing the profile function obeys the ordinary 
differential equation 



Linearizing fl2.6p and requiring a decaying solution at spatial infinity implies an upper limit 
on the frequency uo given by 

l(PU 



2 ^ 2 




9 (2.7) 
2 d/^ /=o 

Integrating (12.61) once, and using the boundary condition at infinity, gives the first order 
equation 

-uo'^f + U. (2.8) 

The boundary condition that the derivative vanishes at x = requires that /(O) = /o > is 
a solution of the equation 

U{h) = u'fl (2.9) 

This equation places a constraint on the type of potential U{f) that allows Q-balls, and 
furthermore, for an allowed potential it determines a lower limit on the frequency uj. This is 
traditionally expressed in the form 

u^>ul=min(^^, (2.10) 

but this assumes that for all / > then U{f) > 0, that is / = is the global minimum of 
the potential. If U{f) < for any / > then clearly once the condition (12. 7p is satisfied 
then the additional condition (12. 9p places no further constraint on the frequency uj beyond 
the original assumption that uj > 0. Obviously, if U{f) < for some / > then it is to be 
expected that some Q-ball solutions will be unstable, since they are nonlinear excitations 
built upon the vacuum = 0, which will be a false vacuum. 

For a potential which is polynomial in the simplest form allowing Q-ball solutions 
is the general 0^ potential 

Uif) = f-f + /3f, (2.11) 

where, without loss of generality, the coefficients in front of the first two terms have been 
scaled to unity using the freedom to rescale the spacetime coordinates and the field. In the 
above potential /5 is a non-negative parameter of the theory. Writing (12. lip in the form 

U{f)=f{l-lfr + {f3-l)f. (2.12) 

makes it clear that if /? > | then / = is the unique global minimum of the potential. If /3 = 
I then there are degenerate vacua at / = and / = ^/2. Q-balls in theories with degenerate 
vacua were first considered in [18]. Finally, if P < ^ then the potential is unbounded from 
below, so / = is a false vacuum. However, even in this case / = is still a local minimum, 
so there will be Q-ball solutions built upon the false vacuum. 

With the general 0^ potential (12. lip the nonlinear field equation (12.50 becomes 

d^dy + 0(1 - 2|0| + 3/?|0n = 0, (2.13) 
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and the profile function equation (I2.6P is 



with solution 



c.V + /(l-2/ + 3/?n, (2.14) 



^'^^ (2.15) 
^1 + VI - 4/5tu'^ cosh(2cj'a;) ' 



where we have defined the complementary frequency uj' = \/l — uP-. 

Using this exact solution the charge and energy can be calculated explicitly to be 
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These explicit expressions reveal that the existence of large Q-balls, that is Q-balls with 
arbitrarily large values of the Noether charge Q, depends crucially upon the parameter /5. 
Large Q-balls exist if and only if /3 > i, which is the requirement that = is the unique 
global minimum of the potential. In this case u;_ < < cu-t-, where the limiting values are 
uj^ = ^Jl — and uj^ = 1. Both the charge Q and energy E are decreasing functions of cu 
in this range, tending to zero as the upper limit uj+ is approached and growing without limit 
as UJ ^ U-. In Figured] the two upper curves display the charge Q (solid curve) and energy 
E (dashed curve) as a function of for /3 = 1/2, in which case u_ = l/V^. 

If /? < I then there is an upper bound on the value of Q for which Q-ball solutions exist, 
as expected in a theory where = is not the unique global minimum of the potential. In 
this case = and = 1, so the lower limit on the frequency is removed, but Q-balls 
are unstable below a critical value u < ujc, as we shall see shortly. In Figured] the two lower 
curves display the charge Q (solid curve) and energy E (dashed curve) as a function of u for 
the value /3 = 1/4. 

The properties of small Q-balls, which have arbitrarily small values of Q, were first 
investigated in [18], and are not very sensitive to the value of /3. Small Q-balls with Q ^ 1 
correspond to u close to the upper limit of the allowed range, a; ~ 10;+ = 1, that is < cu' ^ 1. 
In this limit the above expressions for Q and E can be expanded in powers of u' to give 

Q = 4Lj' + hl6(3-6)u'^ + 0{u'^), E = 4u' + hl6P-8)u''^ + 0{ij'^). (2.18) 

Clearly, the properties of small Q-balls only have a weak dependence upon the value of j3, 
and the leading order behaviour is independent of (3. In fact this leading order behaviour is 
familiar from the well-known limit of the Klein- Gordon equation to the nonlinear Schrodinger 
equation for small amplitude breathers. 

In this paper we shall be particularly concerned with small Q-balls, and the above analysis 
shows that, at least to leading order, the results are universal for all values of f3. We may 
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Figure 1: The charge Q (sohd curves) and energy E (dashed curves) as a function of uj for 
(3 = 1/2 (upper curves) and /3 = 1/4 (lower curves). 



therefore make a convenient choice and set /5 = 0, which we refer to as the truncated modeL 
Taking the hmit /3 = of the above results yields a simplified expression for the profile 
function 

/ = cj'sech(w'a:), (2.19) 

as well as the charge and energy 

Q = Aujuj\ E = (l + 2cj2). (2.20) 

3 

It should be noted that the properties of small Q-balls will be well-approximated by the 
expressions in (12.19^ and (12.20p for the whole class of theories where the potential has an 
expansion in of the form 

Uif)=m'f-g'f + ... (2.21) 

where m and g are positive constants which can be set to unity by suitable rescalings of the 
field and spacetime coordinates. 

Any model with small Q-balls will have the property that Q-balls tend towards fundamen- 
tal particles as the charge tends to zero, with the limit E/Q — > m. The dispersion relation 
for a plane wave solution (p = exp{i{ujt + kx)) of the linearized theory is u"^ = m? + fc^, hence 
uj"^ > and this matches smoothly onto the end of the Q-ball frequency range uj"^ < m? with 
the exponential decay of the Q-ball solution being interpreted as a complex wave number 
k = i\/vp? — uj"^. 
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The condition for a Q-ball to be a classically stable solution is [12] 



^ < 0. (2.22) 

For the truncated model {[3 = 0) then Q = Aoju' and therefore a Q-ball with uj > ^ = ojc'i^ 
stable. It is not surprising that tUc = is the critical value at which instability emerges since 
from fl2.19p this value corresponds to /(O) = which is the turning point of the truncated 
potential U = f"^ — The truncated model has a lower limit on the charge, Q < 2, which 
contrasts with the theory with f3 > where stable Q-balls with arbitrarily large values of 
Q exist. 

2.2 Q-balls in theories with a conformally flat target space 

An interesting generalization of the standard theory (12. ip is to consider a theory in which 
the complex field takes values in a conformally flat target space. Explicitly, the Lagrangian 
density takes the form 

^ = -G^f^^9y-Wm, (2.23) 

where G(|0|) is the inverse conformal factor of the target space metric and W^(|0|) is the 
potential. Obviously the standard Lagrangian (12. ip of the previous section is recovered by 
the choice Cdf/)]) = 1 and the identification iy(|(/)|) = t/(|0|). Again we assume that the 
potential vanishes when the field is zero, that is, W{0) = 0, and the metric is normalized by 
G(0) = 1. 

The conserved Noether charge in this theory is given by 

Q = W (05*0 - 09i0) rfx. (2.24) 

We shall mainly be concerned with the choice G{f) = 1 — /^, which has two main advantages. 
The first is that, as shown below, any potential (which is finite for finite field values) that has 
a mass term yields Q-ball solutions. This contrasts with the standard kinetic term which, as 
discussed in the previous section, requires a rather careful choice of potential so that U{f)/ p 
has a minimum at a positive value of /. Of course, one may take the view that the careful 
choice of U has simply been replaced by a careful choice of G. The second advantage is that 
this choice of G includes an integrable theory, where exact solutions describing multi-Q-ball 
dynamics can be constructed in closed form. 

To begin with, Q-ball solutions will be discussed for arbitrary functions G and W. This 
is most conveniently addressed by making use of the momentum density, as follows. 

The momentum density, p, and the associated current J^, which follow from the La- 
grangian density (I2.23p . are given by 

p = -^{dt<j)dj + 9*09,0), = -^(|9i0|2 + |9,0n + W. (2.25) 
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The conservation law dtP = d^J^ implies that the momentum P is conserved, where 

P= pdx. (2.26) 



A field of the Q-ball form (12. 3p has p = and hence is constant. Therefore, 




where the constant value of is evaluated to be zero by using the boundary condition 
/ — >• as a: — > oo. 

The profile function of the modified theory (12.231) therefore satisfies the first order equa- 
tion 

Note that this equation is identical to the profile function equation (12. 8p after the identifi- 
cation U = GW. 

With the choice G{f) = 1 — then solutions exist for all positive uj satisfying 



2/2 ^d^W 

UJ < ui — 



2 dp 



(2.29) 



since the second constraint (12.91) with U{f) = (1 — P)W{f) is then automatically satisfied 
providing W{1) is finite. 

The simplest theory with G = 1 — p contains only a mass term, which we normalize to 
heW = p. This turns out to be an integrable theory, and we discuss it in detail in the next 
section. Note that in this case U = GW = p — p^ , so the Q-ball solution is identical to that 
of the truncated model with a standard kinetic term. 

Although we have restricted our discussion to theories in (H-l)-dimensions, a similar 
analysis can be performed for the general theory (12.231) in any number of space dimensions. 
The property that a non-standard kinetic term, such as G = 1 — removes the usual 
restrictions on the potential to allow Q-balls is true in any dimension. Non-standard kinetic 
terms are currently popular in cosmological applications and arise in string theory, therefore 
there may be applications of Q-balls in this context to systems not previously considered. 



3 An integrable theory with Q-balls 

The integrable model referred to in the previous section is known as the complex sine- 
Gordon model [3 E] and has a Lagrangian density given by 

-^: = Y^(5.05>-i0r+i0r) , (3.1) 
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where we have chosen to write all terms over a common denominator. The numerator in 
(13.1 p is the Lagrangian density of the truncated model and the two theories share the same 
Q-ball profile function (I2.19p . for the reason described earlier. Later, we shall be interested in 
studying small Q-balls in the complex sine-Gordon model, in which case |0| -C 1. However, it 
is not clear which terms in the Lagrangian can be neglected in this limit, and the fact that the 
profile function agrees with that of the truncated model (for all cj, including small Q-balls) 
suggests that the dominant contribution is obtained by approximating the denominator in 
(13. ip by unity. By comparing solutions we shall see later that this is indeed the case. Note 
that simply keeping all terms to quadratic order in in (13. ip does not provide a good 
description of small Q-balls since such a theory does not allow Q-ball solutions. 
Using the profile function (12.190 the charge and energy are calculated to be 

Q = 4cos"^cj, E = Auj'. (3.2) 

For small Q-balls the leading order behaviour is 

Q = Au' + 0{uj'^), E = 4uj', (3.3) 

which to first order in u' agrees with the result (12.180 for small Q-balls in the standard 
theory with potential (12.111) and any value of (3. 

As uj then naively Q — >■ 27r, however, in this limit /(O) 1 and the integral which 
determines Q is ill-defined and has no unique limit. The integrable model therefore shares 
with the truncated model the property that there is a maximum value of Q, but this is 
realized by two different kinds of behaviour in the two theories. 

The complex sine-Gordon model has been studied by a number of authors, for example [B], 
and several different methods have been used to generate multi-soliton solutions. However, 
it appears that the connection to Q-balls has not been exploited previously, though it has 
been noted that solitons of this theory are examples of Q-balls [15]. Interpretating the 
known multi-soliton solutions in the framework of Q-balls and performing some analysis of 
these solutions will prove instructive in understanding Q-ball dynamics and interactions in 
non-integrable theories. In particular we shall make extensive use of the two-soliton solution 
reproduced below. 

Lorentz boosting the stationary Q-ball solution yields 

(j)k = e'^''w'f,exp(^iwk'yk{t - Vk{x - afc))^sech(^w^7fc(a; - - Vkt)^ , (3.4) 

where 9k is an arbitrary constant phase, Uk is the frequency, cu^ is the complementary fre- 
quency, Vk is the velocity, with associated Lorentz factor = 1/ yl — t^, and is the 
position of the Q-ball at time t = 0. 

Defining the associated complex kink function 

V'fc = -w^tanhf w;^7fc(x - Ofc - Vkt) ) - iwk , (3.5) 
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and the constant 




l+Vk 



(3.6) 



then two Q-ball solutions (j)i,(j)2, with associated kink functions ipi,ip2 can be combined to 
produce the two-sohton solution [6] 



Generically, this solution describes the scattering of two Q-balls, with the traditional inter- 
pretation being that the two Q-balls pass through each other with their individual properties 
preserved and a time advance and phase shift being the only remnant of their interaction. In 
this and later sections we shall analyse this solution in some detail and see that it contains 
much more information than the above simple asymptotic picture might suggest. 

In contrast to most integrable systems there is no problem in setting both speeds vi and 
V2 simultaneously to zero in this two-soliton solution. This yields the solution 



where we have defined X/. = uj[{x — ak) , and without loss of generality we have set = 0. 

The solution (13. 8p describes a breather, that is, a non- stationary solution which is periodic 
in time. The breather frequency is \uj2 — and the solution degenerates to the trivial 
solution = in the limit uJi = uj2- Breathers are usually associated with a soliton anti- 
soliton bound state (for example, in the sine-Gordon model) but here the breather is formed 
from two Q-balls, not a Q-ball anti-Qball configuration, though a breather solution of this 
form can also be obtained by considering negative values of uj2- 

Given that the breather is constructed by setting both velocities to zero in the two- 
soliton solution, then one might expect that this solution describes two stationary Q-balls 
with fixed positions and amplitudes. If the Q-balls are well-separated, that is, the separation 
is much greater than either width, |ai — 02] ^ max(c(j^~^, 0^2"^), then this interpretation is 
approximately true, although for any finite separation the position and amplitude of each 
Q-ball oscillates. The amplitude of the oscillation is a decreasing function of the separation 
between the Q-balls. 

If both constituent Q-balls are placed at the same location (without loss of generality 
this may be taken to be Oi = 02 = 0) then it describes an interesting oscillating periodic 
solution. If uji and 002 are comparable then the oscillation is between a state that resembles 
two resonably separated Q-balls, each with charge 2cos~^co'i -|- 2cos~^to'2, and an excited 
state of a single Q-ball with charge 4cos~^co'i + 4cos~^co'2. Note that because Q-balls are 
non-topological solitons then there is no rigorous definition of the number of Q-balls in 
a given field configuration, so it is perfectly reasonable for a field to have an appropriate 
interpretation in which the number of Q-balls is different at various times. In contrast. 



(^101 - (5202)(^2^1 - ^1^2) - {Sl(j)2 - ^2</'l)(^2'02 - Siljji) 
61 - 5l52(0102 + 0102 + ^^1^2 + V'lV^2) + S2 



(3.7) 



1 — UJ1UJ2 — uj[u2 ( tanhXitanhX2 + cos((u;2 — ^^i)^ + 6'2)sechXisechX2 
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topological solitons have an associated conserved integer that defines the number of solitons 
in any field configuration. Even though a topological multi-soliton solution may resemble a 
single structure, the integer soliton number still allows its identification as a multi-soliton. If 
one of the frequencies is much closer to unity than the other (say a;2 ~ 1) then the oscillation 
resembles a small perturbation of a single Q-ball with frequency uji. We shall see later that 
this breather description of a perturbed Q-ball provides a good explanation of observed 
Q-ball vibrational modes for non-integrable theories. 

In the following section we investigate the force between two well-separated Q-balls in a 
general theory. This will provide some understanding of the fact that the breather solution 
requires the constituent Q-balls to have distinct frequencies. 

4 The force between Q-balls 

Numerical simulations of time-dependent nonlinear field theories in one, two and three 
space dimensions, have revealed universal features regarding the force between two Q-balls 
[31 [2]. The numerical results show that two Q-balls with the same frequency can attract 
or repel depending on whether their internal phases are aligned or anti-aligned, and if the 
frequencies are distinct then a complicated charge exchange process results. In this section 
we provide an analytic calculation of the force between well-separated Q-balls that explains 
this behaviour. 

The force between two well-separated topological solitons is often obtained by using a 
linear addition (or perhaps product) ansatz and calculating the energy of this field as a 
function of soliton separation, in order to derive the interaction potential. However, for 
Q-balls this method can not be applied since the charge Q of such an ansatz will vary with 
the separation, making a comparison of the energy a meaningless result. Explicitly, consider 
the addition ansatz 

= 0, + 02 = e-^Vi + e^-2*+*72 (4.1) 

where /i = fu)i{x + a) and /2 = /i^al^^ ~ with fu){x) the profile function for a Q-ball at the 
origin with frequency uj. Here we choose a to be positive and much greater than the width of 
either Q-ball, so that the field (14.11) is a good approximation to two Q-balls with frequencies 
uoi and uj2 at positions —a and a repsectively, with 9 the relative phase at t = 0. It is a 
simple task to verify that the charge Q of this field is not independent of a, so computing 
the interaction energy as a function of a is not a useful quantity. 

An alternative approach to calculating inter-soliton forces is to identify the force exerted 
on one soliton by the other with the induced rate of change of momentum [13j. We now 
apply this approach to study the force between Q-balls in the general theory (12.231) . which of 
course includes the standard theory (12. ip as a special case. The force can also be calculated 
using the tail interaction methods of [8J. 

Let P[xi,X2\ denote the momentum on the interval [xi,a;2], that is, 

P[xi,X2] = / pdx = - —{dt(f)d:^(j) + dt(f)d^(f))dx. (4.2) 

Jxi J x\ \Jr 
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The force on this interval, X2], is given by the time derivative of the momentum, and the 
conservation law dtp = d^J^ can be used to express this in terms of the pressure difference 
at the end points 



F[xi,X2] 



dP[xi,X2] 

Jt 





X2 




Xl 



-^{\dM' + \dx<P\') + W 



(4.3) 



Xl 



We now apply the ansatz ( 14. ip and note that the interval [xi,X2] of interest will have end 
points which are far from both Q-balls, so that 0i and 02 are small where their evaluations 
are required, allowing the leading order asymptotic force to be calculated by keeping only 
terms up to quadratic order in 0. The functions G(|0|) and iy(|0|) in (14.31) can therefore 
be replaced by G = 1 and W = |0p using the earlier normalizations. Using the fact that a 
single Q-ball has a vanishing current, the leading order contribution to the force is given by 



F[xi,X2] 



- dt(f)idt(j)2 - <9t01<9t02 - <9x01<9x02 - 5x01^x02 + 0102 + 0102 

2 cos (((^2 — 1^1)^ + 6*) 



Xl 



w . dfi df2 



X2 



(4.4) 



Xl 



In order to identify the force on a soliton with the rate of change of momentum on a large 
interval containing this soliton, it is necessary that there is no net flow of momentum in 
an interval which is far from both solitons. For theories where a time-independent soliton 
exists (which is the usual situation for topological solitons) this is automatic, but Q-balls are 
non-topological solitons which are only stationary not static, and therefore this requirment 
is not guaranteed. 

To check the above requirement consider an interval around the origin [—L, L] which is 
far from both solitons, that is, a ^ L. The asymptotic tail of a Q-ball is given by 



f^{x) ~ A^exp{-uj'\x\), 



(4.5) 



where is a positive function of u. Therefore in the interval [—L, L] we can use the tail 
approximat ions 



/i ~ A^^ exp(-u;'i(x + a)), /2 ~ A^^ exp(u;2(a; - a)). 



(4.6) 



Substituting these expressions into (14. 4p yields the result 
F[-L, L] = cos {{uj2 - uji)t + ^)4A^,A^,e-'^KM)(i + uj'^lo'^ - UJ1UJ2) sinh((4 - uj[)L). (4.7) 

This force is identically zero if and only if ui =002, therefore this is a special case in which 
we can calculate the asymptotic force using the above approach. For distinct frequencies the 
result (14. 7p reveals that there is a breather-like motion with a frequency \uj2 — Indeed, 
we have seen that in the integrable model there is an exact breather solution (13. 8p with this 
frequency. 

Restricting to the case of equal frequencies lji = 002 = uj, gives F[—L, L] = and we can 
therefore calculate the force on the soliton at x = a by computing F[L, 00], which will be 
independent of L by the previous calculation. 
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At X = oo there is clearly no contribution to the force because of the exponential decay 
of both /i and /2. At x = L we can again use the asymptotic expressions fl4.6l) but now with 
equal frequencies 



/i ~ exp{-uj'{x + a)), f2 ~ exp{uj'{x - a)). 



Substituting these formulae into the expression (14.41) yields 



F = F[L, oo] = -2cos^ (l-^i^2)/i/2 



dx dx 



-4: cos 9 A^, uo'^e 



x=L 



(4.9) 



which is indeed independent of L, and has the expected exponential fall-off with separation 
2a. 

For two in-phase Q-balls, 6 = and hence F < 0, so there is an attractive force between 
the Q-balls, whereas if the Q-balls are exactly out-of-phase, 6 = n and hence F > and 
there is repulsion. This agrees with the numerical simulations of the nonlinear field theory 
[3|, [2]. If ^ is an integer multiple of tt then, by symmetry, two Q-balls with equal initial 
frequencies will remain equal. However, for other values of 6 there is no symmetry argument 
to keep the initial frequencies equal, and indeed generically a breather-like mode will evolve 
that induces a frequency difference, as discussed in the next section. Therefore the force 
(14. 9 p can only be applied to study the dynamics of Q-balls which are exactly in-phase or 
exactly out-of-phase. 

For both the truncated model and the integrable model, the coefficient Ai^ in the tail 
approximation (14. 5 p is given by A^^ = 2a;'. Hence in these two theories the asymptotic force 
is given by 

F = -lecos^cj'V^'^'^', (4.10) 

which is extremely weak for small Q-balls because of the u"^ factor; recalling that the Q-ball 
mass approaches iu'. 

The dynamics generated by the asymptotic force (I4.10p is obtained by solving the equa- 
tion Md = F, where the mass M is the energy of a stationay Q-ball with frequency u. 
Restricting to the truncated model then the mass is given by the second expression in (I2.20p 
and the equation of motion for two in-phase Q-balls becomes 



l + 2cj2 



(4.11) 



Taking initially stationary Q-balls, a(0) = 0, with a(0) = Oq, the solution to this equation is 



1 , / /2V3cj'2e-''«^'t 
a = an^ log cos , — 



(4.12) 



As an example, for the values uj = 0.98 and oq = 25, the position (I4.12p is plotted as the 
dashed curve in Figure [21 for times until the position vanishes. The solid curve on the same 
plot is the result of a numerical simulation of the full nonlinear field theory with the initial 
conditions formed using the addition ansatz (14. ip with the same parameter values. In the 
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Figure 2: The Q-ball position as a function of time, for an in-phase two-Q-ball system with 
equal frequencies uj = 0.98 and initial positions ±25. The solid curve is the result of a 
numerical simulation and the dashed curve is the approximation based on the asymptotic 
force. 



field theory simulation the position of each Q-ball is computed as the point at which \(f)\ is 
maximal. Once the two Q-balls are sufficiently close together (a ~ 6.5 in this example) the 
field 101 no longer has two distinct maxima, but rather has a single maximum at the origin. 
At this point the two Q-balls lose their individual identities and it does not make sense 
to assign two Q-ball positions, which is why the solid curve terminates in Figure [2l The 
evolution of the system beyond this point follows a damped breather motion, as discussed in 
the following section. There is excellent agreement between the two curves in Figure [21 even 
down to the smallest allowed separation. This demonstrates the accuracy of the asymptotic 
force and reveals that it provides a good description well beyond the expected region of 
validity. 

Returning to the case of distinct frequencies, uji uj2, then although the force on a small 
interval far from both solitons (14.71) is non-zero, it does vanish if averaged over the time 
period T = 27t/\uj2 — uJi\. It would therefore make sense to compute the time-averaged force 
on a Q-ball, but the expected leading order contribution also averages to zero. We therefore 
see that two Q-balls with distinct frequencies perform a breather-like motion where the 
expected leading order contribution to the force averages to zero. In fact, as we have seen, 
the integrable model has an exact breather solution, hence in this theory the average force 
vanishes to all orders in the separation. For a generic theory one expects a higher order 
contribution to be non-zero, which implies a very weak force between Q-balls. 

Note that the fact that the force does not vanish in the integrable model if the two 
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frequencies are equal explains why this limit degenerates in the exact breather solution. 

Similar calculations to the ones performed in this section, to determine the rate of change 
of momentum, reveal analogous results for the rate of change of charge. The breather-like 
motion is therefore also associated with a quasi-periodic flow of charge and relates to the 
charge exchange observed in numerical simulations [3j. This is discussed further in the 
following section. 

5 Breathers and multi-Q-ball dynamics 

The exact breather solution (13. 8p of the integrable model can be used as an initial condition 
in the truncated model, to examine the differences between integrable and non-integrable 
evolution. The expectation is that the breather motion is damped, as the non-integrable 
model radiates energy into other modes. A typical simulation is presented in Figure [3] with 
the initial conditions obtained from the integrable solution with parameters ui = 0.99, uj2 = 
0.98, 02 = TT, ai = 02 = 0, corresponding to a period of T ^ 628. The graph displays 
the distance d of either Q-ball from the origin, where the Q-ball position is calculated as 
the location at which |0| is maximal. The inset graph shows the evolution for the first few 
oscillations (0 < t < 4000) revealing a slow decay of the amplitude, indicating that the 
maximal separation of the two Q-balls slightly decreases with each cycle. Note that each 
oscillation contains a time interval where the separation identically vanishes. During such 
an interval the maximal value of |0| performs a single oscillation, but the maximal value 
remains located at the origin, that is, the configuration contains a single peak rather than 
two peaks during this time interval. 

The main graph is a plot for a much longer time period, < t < 35000, and confirms that 
many oscillations are performed before the system becomes a single Q-ball configuration, in 
the sense that there is a single peak for |0| for all time. In fact, a plot of the maximal 
value of as a function of time shows that the configuration continues to oscillate, with a 
reasonably large variation in |0| at the origin, for a length of time which is far beyond that 
shown in Figure [31 This part of the evolution also has an approximate breather description, 
but with one of the frequencies much closer to unity than the other. Recall that in this 
parameter regime the exact breather solution of the integrable theory also has the property 
that it remains a single peak throughout the entire period. This will be analysed further 
when we discuss vibrational modes of a single Q-ball. 

The dynamics described above may therefore be viewed from two different perspectives. 
The first description is an oscillating two Q-ball system which dissipates energy until it 
becomes an oscillating single Q-ball and eventually a stationary single Q-ball. The second 
description is as a slow drift through a family of breather solutions (labelled by the pa- 
rameters Ui and U2) in which one frequency tends towards unity while the other frequency 
decreases to approximately conserve charge. Taking the limit co'2 — 1 in the breather solution 
(13. 8p yields the stationary single Q-ball solution with frequency ui, so both descriptions are 
consistent. 

The single stationary Q-ball which results from the above evolution has a charge which is a 
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Figure 3: The distance d from the origin of the maximum of |0|, for times < t < 35000. The 
simulation is performed in the truncated model. The initial conditions (field and its time 
derivative at t = 0) are taken from the exact breather solution (13.81) of the integrable model, 
with uji = 0.99, UJ2 = 0.98, 02 = vr, oi = 02 = 0. The inset is a blow-up of this plot for early 
times < t < 4000. There are many oscillations with a slowly decreasing amplitude. 



few percent less than the initial configuration. The charge appears to be lost gradually, rather 
than in sudden drops, but this is difficult to accurately measure numerically, particularly 
because it is a small effect. Moreover, we are unable to present convincing evidence to 
discriminate between the charge and energy being lost through radiation or through the 
emission of tiny Q-balls. We suspect that the latter is correct, but because Q-balls tend 
towards elementary particles as the charge tends to zero, this makes it tricky to distinguish 
the two in a noisy background. 

In the numerical simulations described above the spatial domain has a finite range. If 
Dirichlet boundary conditions are applied than charge and energy are conserved throughout 
the simulation to a very accurate precision. However, with such boundary conditions energy 
and charge are observed to propagate to the boundary where they are refiected back into 
the simulation domain. This is a standard phenomenon with Dirichlet boundary conditions 
and does not accurately represent the true dynamics on an infinite domain. To minimize the 
refiection produced by the boundary we apply absorbing boundary conditions. This consists 
of a small region close to the boundary where a damping term is introduced to dissipate 
most of the energy and charge before it is refiected back into the bulk of the simulation 
region. This provides a better approximation to the true system on an infinite domain, and 
in particular allows charge and energy to decrease, as it should for the theory defined on the 
whole line when energy and charge are only computed in a finite interval. 

Similar results to the one described above are obtained for other parameter values. As 
the charge of the configuration increases the damping becomes more significant and less 
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Figure 4: Plots of at increasing times in a two Q-ball scattering process with frequencies 
uji = 0.98, UJ2 = 0.99, initial positions ai = —02 = —25, and speeds vi = —V2 = 0.1. Solid 
curves are in the truncated model and dashed curves (often hidden by the solid curves) are 
in the integrable model. 



oscillations are required until a single peak forms. This is to be expected, since the integrable 
model is a better approximation to the non-integrable theory for smaller charge. Note that 
the range of breathers that can be studied in the truncated model is restricted, since stability 
requires that the field must satisfy the constraint |0| < for all points in space at all times. 
During a breather cycle the maximal value of |0| can be substantially larger than it is for 
either two distinct Q-balls or a single Q-ball, even when the configuration has the same 
charge as the breather. 

The evolution of two equal frequency in-phase Q-balls, discussed in the previous section, 
also yields a damped breather motion, once the two Q-balls merge to form a single peak 
for the first time. As the breather solution (13.81) degenerates for equal frequencies, the 
appropriate breather parameters, uji and uj2, to compare with (at least when the total charge 
is small) are obtained by matching the total charge and energy of the initial configuration 
to that of the breather. 

Recall that the two-soliton solution (13.71) describes the elastic scattering of two Q-balls 
which emerge with their charges and speeds preserved. We therefore expect the scattering 
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of sufficiently small Q-balls to be close to this integrable behaviour for all systems in our 
generic class of theories. 

As an example, the solid curves in Figure H] are plots of |0p at increasing times in a two 
Q-ball scattering event in the truncated model. The initial condition is generated using the 
simple addition ansatz with Q-balls of frequency ui = 0.98 and uj2 = 0.99 and respective 
positions ai = —02 = —25. Both Q-balls are initially moving towards each other with equal 
speeds, vi = —V2 = 0.1, and are in-phase (6' = 0) at t = 0. For comparison, the same 
scattering event in the integrable theory is shown as the dashed curves in Figure H] (though 
the two evolutions are so close that it is difficult to distinguish the sohd and dashed curves). 
Note that the exact solution (13.71) could have been used to plot the dashed curves in Figure HJ 
but as a check on the numerical code and initial conditions, the integrable theory was also 
solved numerically - reproducing the elastic scattering result to a very high precision. 

The simple asymptotic picture of trivial scattering in the integrable theory hides the 
complicated interaction which takes place when the Q-balls coalesce. During this part of 
the scattering process the configuration is highly distorted and the evolution resembles a 
breather-like motion. In the integrable theory the two Q-balls emerge at exactly the point 
in the breather cycle at which the individual charges are identical to the incoming charges. 
In the non-integrable theory the breather is slightly distorted and the Q-balls emerge at a 
point in the cycle at which the charges differ slightly from the initial values. With increasing 
charge the breather cycle and the point of separation diverge between the integrable and 
non-integrable models. This leads to a significant difference between initial and final charges 
in non-integrable models, and results in the charge exchange phenomenon seen in earlier 
numerical simulations [3]. 

6 Perturbing a Q-ball 

Perturbations of a single stable stationary Q-ball lead to extremely long-lived oscillations. 
In this section we discuss such oscillations and provide formulae for the frequencies of the 
vibrational modes excited. 

Given a Q-ball solution (12. Sp of a standard theory (12. ip . we wish to consider perturbations 
which preserve the total charge (12. 2p . and therefore allow the perturbed solution to eventually 
return to the initial unperturbed stationary Q-ball. A simple example is a modified squashing 
perturbation, which uses an initial condition taken from the field configuration 

^ = v^e*-*/(Ax). (6.1) 

If A = 1 then (16.11) is the stationary Q-ball solution, but for other positive values of A it 
describes an initial squashing or stretching of the Q-ball. It is easily checked that the charge 
(12. 2p of this initial field configuration is independent of A, and hence this is an example of a 
charge-preserving perturbation. 

Figure [5] displays part of the evolution for an initial squashing of the above form with 
A = 1.05 and Q-ball frequency w = ^ in the truncated model. The plot shows the amplitude. 
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Figure 5: The amplitude, that is |0| at x = 0, as a function of time t G [700,800], for a 
perturbed Q-ball with frequency u = ^ and squashing parameter A = 1.05. 
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Figure 6: The power spectrum of the oscillation displayed in Figure [5l 
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that is 101 at X = 0, for the time period 700 < t < 800. This demostrates the almost quasi- 
periodic nature of the oscillation and the very slow decay of the amplitude. 

Figure E] is a plot of the power spectrum associated with the oscillation presented in 
Figure O It is clear that there are two dominant vibrational modes which are excited, 
with frequencies Qb = 0.132 and Qp = 1.516 in this example. Other charge-preserving 
perturbations yield similar results, as do scattering processess in which a Q-ball is excited 
through interactions with other Q-balls. In each case two dominant vibrational modes are 
excited and the aim in this section is to understand these modes and obtain formulae for the 
frequncies Qb and flp a.s a. function of the initial Q-ball frequency uj. 

Recall that the exact breather solution (13.81) of the integrable theory has a frequency 
\u!2 — where Ui and 002 are the frequencies of the two constituent Q-balls. Since the 
energy of a Q-ball tends to zero as the frequency tends to unity, any arbitrary perturbation 
of a single Q-ball contains enough energy to create an additional Q-ball. Assuming the 
perturbation is small then any created Q-ball must have a frequency close to unity. Taking 
ui = u, the frequency of the initial unperturbed Q-ball, and approximating the frequency of 
the created Q-ball as ~ 1 then we may interpret the perturbed Q-ball as a breather with 
frequency Qb = ^2 ~ ^1 ^ ^ — ^- the numerical simulation presented above the values 
give 1 - = 1 - ^ = 0.134 « 0.132 = ^b- 

In summary, we see that the dominant mode excited in the perturbation of a single Q- 
ball is a breather mode, with frequency f^e ~ 1 — t^- Other simulations of perturbed Q-balls 
yield similar results and are consistent with this formula and its interpretation. In particular 
this interpretation explains why the mode is so long-lived, as it is close to a solution of the 
integrable theory which is strictly periodic. Note that if the perturbation is large then the 
energy imparted to the created Q-ball can be significant, so the approximation that 002 ^ 1 
may no longer be valid. In that case the frequencies uji and uj2 should be calculated by 
matching the energy and charge of the perturbed Q-ball to the breather with parameters 
uJi and UJ2- The result will be that 002 is slightly less than 1 and uti is slightly greater than 
a;, in order to preserve total charge. Thus the difference Qb will be slightly reduced in 
comparison to the simple approximation Qb ~ 1 — w. This has also been confirmed in 
numerical simulations. 

The interpretation of a perturbed Q-ball as a breather is similar to a recent observation 
by Kalberman [9] on the oscillation of a perturbed kink, called a wobble, in the sine-Gordon 
model. In that case a perturbed kink can be described by an exact solution consisting of both 
a kink and a breather component. The difference for Q-balls is that the original soliton is 
non-topological, so a closer analogy in the sine-Gordon theory would be the perturbation of a 
breather solution described by a two-breather solution. The generic phenomenon is that any 
theory containing a non-topological soliton with arbitrarily small energy will produce such a 
solution as a dominant low energy vibrational mode when any other soliton (topological or 
non-topological) is perturbed. If the model is integrable then an exactly periodic oscillation 
will result, but if the theory is close to an integrable theory then the perturbation will yield 
a long-lived oscillation with a slowly decreasing amplitude. 

Now we turn to the second mode, with frequency Qp, observed in the perturbation of a 
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Q-ball. We shall show that this mode corresponds to an oscillation in the internal frequency 
of the Q-ball by using the moduli space approximation (for a review see [H]) to determine 
the dynamics of the internal frequency by reducing the field theory to a finite dimensional 
system. 

Consider a field configuration of the form 

= a(t)e^^W/(x), (6.2) 

where a and Q are time-dependent real moduli and j[x) is the profile function of a Q-ball 
with frequency lo. Substituting the field (16. 2p into the Lagrangian of the truncated model 
and performing the spatial integration yields the effective Lagrangian 

L = (d^ + ^2^2)2^/ _ ^2^J ^ 4 4^,3 _ 2 2^,3^ ^^^^^ 

o o 

The equations of motion which follow from this Lagrangian are 

a^ + 2d^ = = -^(a^^) (6.4) 

OL at 

a-a9^ + a{l + ^u'^ - ^aV^) = 0. (6.5) 

Equation (16.41) represents the conservation of Noether charge associated with the field (16. 2p . 
This reveals that at least two degrees of freedom must be included in the moduli space ap- 
proximation. If the motion was reduced to only the internal phase 6 then charge conservation 
would force its velocity to be constant, and therefore not allow an oscillation of the frequency 
9. Including the degree of freedom a allows a change in the amplitude to compensate for a 
variation in the frequency. 

The stationary Q-ball solution corresponds to 6 = uj and a = 1, which is easily verified 
to be a solution of the equations (16. 4p and (16.51) . To study small oscillations around this 
solution we set 6 = ut + rj and a = 1 + e and linearize the equations of motion in r] and e to 
obtain 

g 

ii + 2uje = 0, and e - 2ujfi - -u'^e = 0. (6.6) 

3 

Integrating the first equation gives fj = —2uje + c, where c is a constant determined by the 
initial perturbation. Substituting this into the second equation yields 

4 

e + -(5cj2-2)e-2cuc = 0. (6.7) 

Hence the oscillation frequency is given by 

2 



A = -^V5cj2 _ 2, (6.8) 
v3 

which is automatically real since stable Q-balls only exist in the truncated model for a;^ > |. 
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The earlier numerical simulation dealt with a Q-ball with frequency = which when 
substituted into the formula (16. 8p produces 

A = = 1.528 ^ 1.516 = (6.9) 
V ^ 

Perturbations of Q-balls for a range of values of uj confirm that the expression (16.81) is a good 
approximation to the frequency VLp oi the second observed vibrational mode. Furthermore, 
an initial perturbation of the form (16. 2p can be used as a charge-preserving perturbation in 
the full field theory and a resulting power spectrum analysis shows that more of the energy 
is imparted into the VLp mode than for other forms of perturbations, such as (16. ip . 

The results of this section have demonstrated that a perturbation of a stable Q-ball 
excites two main vibrational modes; the first a breather mode associated with the creation 
of a small Q-ball, and the second an oscillation of the internal frequency and amplitude of 
the Q-ball. Explicit formulae for the frequencies of these vibrations have been derived for 
the truncated model and verified to be in good agreement with numerical results. The same 
method can be used to derive an analogue of the formula (16.81) for the general 0^ model with 
potential (12. lip though the resulting expression is more complicated. 



7 Duality 

It is known that the integrable complex sine-Gordon model (13. ip has a dual description 
in which Q-balls of the original theory map to kinks of the dual theory. A reformulation of 
the theory in terms of a gauged Wess-Zumino-Witten model reveals that this is a T-duality 
relation [17], and this has been extended to a range of integrable theories [15]. In this section 
we show that a similar duality transformation exists for any (l+l)-dimensional theory with 
Q-balls and does not rely on the integrability of the theory. 

We shall define a dual theory, with field ip, associated with the general theory (I2.23P with 
field (f). The first part of the construction requires a relation between |'?/'| and |0| given by 
solving the equation 



dW G{\ct>\y 
with the condition that = 1 when |0| = 0. Explicitly, 



(7.1) 



Given this relation we define the dual functions Gd'?/'!) and W^(|-?/'|) via the formulae 

G(iv^i) = ^^^, wm) = wm- (7.3) 

The Lagrangian density of the dual theory is taken to be 

C = -^^d,,l^d^i,-Wm). (7.4) 
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As an example, the complex sine-Gordon model has G(|0|) = 1 — |0p and = In 

this case 

\^\' = ew -T^ = l-|0r, (7.5) 



1-/ 
therefore 

and 

wm) = wm = \<p\' = i-\^p\'. (7.7) 

The dual theory for the complex sine-Gordon model is then recovered as 



1 - , , , 

The complex kink field (13.51) . used in constructing multi-soliton solutions of the complex 
sine-Gordon model, is a solution of the dual theory (17.81) . and we shall discuss this aspect 
more generally below. 

In the original theory the vacuum field upon which Q-balls are built is (f) = and via 
(17. 2p this maps to I?/'! = 1 in the dual theory. Therefore solutions of the dual theory are 
complex kinks which connect two different points on the unit circle as x — >■ ±oo, with the 
phase difference in ip being a real- valued topological charge. We shall denote this topological 
charge by Q and normalize it as 

— ^^^^ — —dx. (7.9) 

Note that the value of ip at spatial infinity is fixed for dynamical reasons since as ^ oo 
then I'i/'l — >■ 1 and therefore by (17.31) G ^ 0, freezing the dynamics of the field at infinity. 

In the following we describe how solutions of the dual theory are obtained from those of 
the original theory. The relation (17.21) is used to obtain lipl given |0| but we still need to 
specify how to obtain the phase of ip. 

The charge density q of the conserved Noether charge Q (12.241) in the original theory is 

(l=77^,i<P9t^-^9t<l)), (7.10) 
G{\(f)\) 

and the associated conservation law is d^q = d^J^ where 

^^ = ^^(05x0-05x0). (7.11) 

This allows the definition of a field G through the relations 

dxQ = \q, and 9*6 = (7.12) 
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We now prove that 

^ = e-^®|^|, (7.13) 

is a solution of the dual theory. 

Combining (17.131) with 07.121) gives 

ijjd^ijj - ipd^j) _ (pdjcp - (j)dt(p jfdttp - ipdtip _ (pd^cj) - (pd^cp 

i^p - cm) ' i^p - GM) ■ > 

Calculating (9^|V^P and using (17.11) to express this in terms of yields 



(7.15) 



Equations (17.141) and (17.151) are most conveniently expressed using light cone coordinates 
M = |(t + x) and V = ^(t — x) to give 

In light cone coordinates the field equation which follows from the dual Lagrangian density 

dzaD is 

Using the expressions in ( 17. 161) together with ( 17.11) and the definitions ( 17.31) for G and W it 
can be shown that 

G(l^l) d\ip\^ d\ilj\^ 

a„.,^a„,a>X^ + ,G(l««|. (7.19) 



The vanishing of the final expression in the above is equivalent to the field equation of the 
original theory in light cone coordinates, hence we have proved that the definition ( 17.131) 
with (17. 2p produces solutions of the dual theory from solutions of the original theory. 

As —0 is the phase of ip the first relation in (17.121) shows that the topological charge 
( 17.9P in the dual theory is equal to the Noether charge in the original theory since 

d,edx= / qdx = Q. (7.20) 

A stationary Q-ball in the original theory has = 0, therefore by the second relation in 
( 17. 12p is independent of time. Therefore we have shown that in any theory a stationary 
Q-ball solution has a description in a dual theory as a static kink, with an interchange of 
Noether and topological charges. It can also been shown that the energy and momentum is 
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preserved under duality and that the Lagrangian densities of the original and dual theories 
are equal up to a total derivative. 

The dual theory (17.81) for the complex sine-Gordon model has kink solutions with a 
continuous range of energies, because the vacuum points at spatial infinity are any two 
different points on the unit circle. It is interesting that embedding these kinks as domain 
walls in a higher dimensional theory therefore yields domain walls with a continuous range 
of tensions and allows the construction of BPS junctions with arbitrary angles [TB] . 

It remains to be seen whether the dual description of Q-balls can be exploited to fur- 
ther investigate the dynamics and interactions of Q-balls, but it may prove useful as static 
kinks are generally easier to deal with than Q-balls, because the complication due to time 
dependent phases does not arise. 

It is interesting to note that there is a self-dual theory defined by 



and is another known example of an integrable theory, generally referred to as the complex 
sine-Gordon II model [7]. This model has degenerate vacua at |0| = and |0| = 1. From our 
earlier analysis the integrable model (I7.2ip shares the same Q-ball solutions as the standard 
theory (^J^ with the potential f/(|0|) = \(f)\^{l - The integrable model (TTm could 

therefore be used to study this theory, in the same way that the integrable model (13. ip was 
used to study the truncated theory. However, as we have already demonstrated, the proper- 
ties of small Q-balls in both integrable theories will converge to the same limit as the charge 
decreases (once the appropriate rescalings have been applied to match normalizations). 

8 Conclusion 

In this paper we have performed some of the first analytic studies of Q-ball dynamics and 
interactions. We have shown how an integrable theory can be used to study small Q-balls in 
non-integrable theories and explain some of the phenomena observed in previous numerical 
investigations. Although most of our analysis has been applied to small Q-balls, and has 
been restricted to (l+l)-dimensions, there appears to be much in common with numerical 
simulations in higher dimensions and for larger Q-balls [3l[2]. This suggests that the results 
should be of some relevance to more general issues regarding Q-balls. One obvious difference 
is that in (3+l)-dimensions a theory with a standard kinetic term and a potential which is 
polynomial in |0p does not allow small Q-balls, although a potential which is polynomial 
in 101 does have small Q-balls [10]. This difference can be understood by examining the 
conditions required for the existence of small Q-balls for a general potential in an arbitrary 
number of space dimensions [L9l . 

Finally, the integrable theory also provides an effective approximation to small Q-ball 
anti-Qball dynamics. As the integrable theory has an exact periodic solution of this type 
then, as expected, in non-integrable theories the behaviour is a perturbation of this periodic 
case. The Q-ball and anti-Qball oscillate together for a number of cycles and the annihilation 




(7.21) 
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process is slow, with only a small amount of charge lost in each cycle. The fact that Q- 
ball anti-Qball annihilation is a slow and inefficient process could have repercussions in the 
cosmological context of Q-ball formation in phase transitions. 
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